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Abstract 

The direct treatment of the Faddeev equation for the three-boson system in 3 dimensions is generahzed 
to nucleons. The one Faddeev equation for identical bosons is replaced by a strictly finite set of coupled 
equations for scalar functions which depend only on 3 variables. The spin-momentum dependence occurring 
as scalar products in 2N and 3N forces accompanied by scalar functions is supplemented by a corresponding 
expansion of the Faddeev amplitudes. After removing the spin degrees of freedom by suitable operations 
only scalar expressions depending on momenta remain. The corresponding steps are performed for the 
deuteron leading to two coupled equations. 

PACS numbers: 21.45.-v, 21.30.-x, 21.45. Be 
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I. INTRODUCTION 

With the advent of chiral three-nucleon (3N) forces at N^LO [l| which have a rich spin and mo- 
mentum structure, the calculation of the 3N system based on standard partial wave representations 
is rather tedious. Moreover, for higher angular momenta it is a challenging task to reliably control 
the numerical accuracy Q]. For a system of three bosons a numerical treatment working directly 
with momentum vectors for bound and scattering states including model three-boson forces turned 
out to be very feasible 

mm 

. Therefore we suggest to also formulate the Faddeev equation 
for the 3N system directly with vector variables. Our proposed formulation is an alternative to 
Ref. [7], where spin and iso-spin couplings are explicitly carried out. The key ingredient thereby 
is an operator representation of the ^He wave function given long time ago by and later on 
modified and strictly re-derived from the complete partial wave representation [9|. Thus, the idea 
is to represent 2N t-operators, the 3N forces and the Faddeev components of the 3N wave function 
as products of scalar functions with scalar products of spin operators and momentum vectors. Fur- 
thermore, by suitable operations one removes the spin operators leading to relations among scalar 
functions of momentum vectors only. 

In Section |TI] we consider the deuteron as warm-up exercise, whereas the central topic, the 
treatment of the 3N Faddeev equation, is given in Section lllli We end with a brief summary and 
outlook in Section llVi 



Mp) + ( '7(2) • p '?(3) -p-^p^ ) Mp) 



II. THE DEUTERON 

The deuteron momentum wave function has the well known operator form (see for instance [lo| ) 

= ^4>k{p)bk{ff(2),^i3),P)\'^1TT'd), (1) 
fc=l 

where {Imd) describes the state in which two spin 1/2 states are coupled to total spin 1 and 
magnetic quantum number m^. The two scalar functions <j)i{p) and </>2(p) are the s- and d-wave 
components of the deuteron wave function. We label the two nucleons as particles 2 and 3, and 
the operator (T(j) (i=2,3), represents the corresponding spin operators. 

As is well known, rotational-, parity-, and time reversal-invariance restrict any NN potential V 
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to be formed with six terms ll|], so that the most general NN force has the form 



(2) 



with Vj{p,f/) being scalar functions. The spin-momentum scalar operators Wj{^(2)i^{'i)^Pii^) ^'I's 
given by 

Wl = 1 

W2 = ('7(2) + o'(3)) -{pxp) 

W3 = di2) -{pxp) (T(3) -{PXP) 

W4 = CT(2) ■{p+p') (T(3) ■{P + P') 

W5 = S(2) -ip-p) O'CS) -{P-P) 

W6 = ^(2) -{P-^) ^{3) -{P + P) 

+ ^(2) ■{P+P') '?(3) ■{p-p'), (3) 

and represent the maximum possible set of invariant operators. Using these operators the 
Schrodinger equation in integral form reads 

1 

'"(2) • P 0-(3) 'P~^i 



h{p) + [^(2)-P^(3)-p-nP ]Mp) 

6 

dV Yl ^) (^(2) ' '^(3) ' ^' 
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We remove the spin dependence by projecting from the left with {lrnd\bk{(T(2)j'^{3)iP) summing 
over md- This leads to 
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J2 {'^^d\bk{ff(2),^{3),P) X] (f'k'{p)bk'{^{2)^(3)P)\'^rnd) 



md=-l 



Ed 



k'=l 
1 ^ 6 



m md=— 1 j=l 



X 0fe"(p') bk"{d(2),^(3),P)\'^T^d) ■ 
k"=l 

If we define the scalar objects 

1 

Akk'ip) = Y (lm'd|&fe('?(2),0'(3),p)6fc'('^(2),<?(3);P)|l"i<i) 
ma=-l 
1 

Bkjk"{P,0) = X] '\'^^d\bk{S(2),S{3),P)Wj{a(2),d{3),P,P) 



(5) 



(6) 



^'fe"(^(2),^(3),^)|l"^d), (7) 

which can be calculated once for all we obtain 

2 ^ „ 6 2 

X] ^kk'ip)(Pk'{p) = p2 / d^p'^Vj{p,f) X Bkjk"{p,^)<t>k"{p') ■ (8) 

k'=l m j=l k"=l 

This is a set of two coupled equations for tlic s- and d-wave components (f>i{p) and (p2{p)- The 
summation over guarantees the scalar nature of A^k'ip) and B^jk/'ip,!^)- Expressions for 
Akk'ip) and Bkjk"{p,i^) are given in the Appendix. 

The azimuthal angle can be trivially integrated out, leading to the final form of the deuteron 
equation 

2 

X] ^kk'{p)(t>k'{p) = 
k'=l 

^ '•' ^ '' ^ ^ (9) 

m k 

where x = p-0. 
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III. THE 3N BOUND STATE 

The Faddeev equation for the 3N bound state reads 

= GotP^ + (1 + Got)GoV^^^ (1 + P)V' , (10) 

where Go is the free propagator, t the NN t-operator for nuclcons 2 and 3, the permutation operator 
P = P12-P23 + -Pi3-f23- The term V^^^ is that part of the 3N force, which is symmetrical under 
exchange of nucleons 2 and 3. The choice of the subsystem pair 2, 3 is of course arbitrary. The 
Faddeev component is given by ip, and the total 3N state is then = (1 + P)'tp- 
We introduce the three possible 3N isospin states 

l7o) = I (ol) h 



^ 2J 2' 

M = I (4) i> 

M = I (li) 2), (11) 

with the notation that the 2N subsystem isospin is either or 1 and coupled with the isospin of 
the third nucleon to the total isospins T=l/2 or 3/2. 
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We expand the Faddeev component ip, the 2N t-matrix t, and V^^^ as 



(12) 



t = El^)*7,7'(7'l (13) 



7:7 



1^^^^ = El^)^7,7'<^'l- (14) 

7:7' 



Then projecting Eq. (jlOp from the left, we obtain 

= Go 5^ t^,y(7'|P|7")V'. 



+ Yl (7l(l + Got)|7')Go^-;/i.(7"l(l + ^')l7"V7"' • (15) 



We neglect the change of the total 2N isospin and the change of the total 3N isospin in the 3NF 
and define 

*7,7' = ^tt'ttTT' (16) 

yi^' - Stt'V,% . (17) 

Furthermore [l^ . 

(7|P|y) = 5rr,p,i,r = 6TT'Fu'T{P!^Pir + (-)*+*>iTAT) , (18) 

where -Ftf'T are geometrical factors and PfJ^ acts only in spin and momentum space. Inserting all 
that into Eq. (fT5]) yields 

T' t' 

+ (<5tt' + GottTT') Go Vi^.^' 

T' t' 

5^(5t't" + Ft,t"T'{PtTP2T + (-)*'+*'>iT^'2T))V't"r' . (19) 
t" 

The charge independence and charge symmetry breaking in the 2N force leads to the coupling of 
total isospin T=l/2 and 3/2. 

In the space spanned by the two Jacobi momenta p and q we encounter the off-shell NN t-matrix 



\ttTT¥i^) = S{q- q')ttTT'{pJ, E,) (20) 
with Eq = E — and the matrices for the permutations [l^ 

\plTP^TWi^) = s{p-n<i,Q'mp-^'iQ,Q))P!2Pi3 

]PtTPiTW^) = 5{p + m q)W + ^'{q, Q))P!3Pi3 (21) 



with 

1 



^%^) = l^q + q (22) 



^'{q,^) = q+^^ . (23) 
Consequently with iptrip, q) = (wIV'tr) we get 

m 4m" 

T' t' 

+ {-f+''ttTT'{p, -aq, ^).Eg) p^pk Yl ^tt'T' #T'(-^'(?, q'),^)) 

t' 

E m Ami t' 



111 ^'in -•■ 1 

J2 [ d'p"dh"v}^},,{j/,qJ',f)^t'T'{j/',f) 
t' 

+ 3 j d\'d\" 

Y (ni'T(P' q"), q") P^Pk E ^t't"T ^f'Ti^'io', q"), q") 

t' t" 
+ i~ f+'"vSl{p, q, -7Ti<f, q"),q") Pf^Pi, ^ F,,^,,^ A"T{-^'{'f, g"),g")) 

t" 

1 f \ 

m 4m ^ T' m 4my 

/ d^q'dV Yl {vSU^^Q^nQ',f),Q'')P!2PisY^t't''T'^^^^^^ 
'' t' t" 

+ (-)*'+*'VS'(p',9,-^(g',9"),r') PIA 

Y Pt'f'T' i^t"T"{-^'{^, g")> ^')) ■ (24) 
t" 

The antisymmetry of the Faddeev component under exchange of particles 2 and 3 imposes the 
condition 

{-fPiT^tT = ^tT (25) 

which in momentum space reads 

{-fPk^tT{-p. q) = ^tT (P, q) . (26) 
Thus V'oT (V'it) is even (odd) under exchange of particles 2 and 3 in spin and momenta. 
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This relation can be used to simplify Eq. (j24p . We regard 

P^Pki^txi-p, q) = Pi-sPuPisPis^tri-p, q) = PkPl2Pk{-fi'tT{p^ 0) (27) 
which turns Eq. ([21]) into 

i'tT{p,q) = 

1 



p. 3 2' 



23; 



t' 

+ 3 , / d^p'd\'Y.V^^)r{p,q,j^ i;t'T{0 ,(^) 

E — Trr — -TTrrQ J xi 



+ 



E- 


m 4m " 




1 


E- 


m 4m " 


E 




t' " 






1 



1 



m 4m ^ 



m 4m ^ 

Tf (g', r), r) + (-)*vi'^(p, -^g, r),g")p 

t' 

PliPk ^t't"T i^f'Ti^'iq', q'), q) 
t" 

H 9 / C^V ttTT' iP, 0, En) 5 

-C' — t;7 — i:;;:? t' — — — tztQ 



t' 

PuPis Yl ' • (28) 

t" 

We now come to the decisive points. In Ref. Q| we derived the operator form of the 3N bound state 
wave function based on the complete set of partial wave states. Exactly the same arguments can 
be applied to the Faddeev components iptrip, q)- Thus, in the notation of Ref. [sl one can write 

8 8 

i^tT{p,q) = Y't>^kp^q) 0,\xn = Y4>fkp^q)\x^), (29) 



i=l i=l 



where |x™) = 1(02)2"^) ^ specific state in which the three spin-1/2 states are coupled to the 
total angular momentum quantum numbers of the 3N bound state. The functions (p[j^{p,q) are 
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scalar functions and the operators Oj are given as 



d 


= 1 












" (1) vi^ ^ y/ 




O4 


= a{23) ■ p X q 




O5 


= (t(23) • q (7(1) 




Oe 


= (t(23) •p<T(i) 


• 


O7 


= (t(23) 




Os 


= (t(23) • g (7(1) 





Here (7(23) = i(a(2) - (7(3)). 

The second expansion in Eq. (f29l) is based on 



X2 = — (7(23) -(7(1)1 



XI = \xn 

1 

1 



X4 = -^(^^(7(23) X (pxq)- (7(1) X (7(23) • {p x q)j\x''- 



X5 = 7 ('^(23) - ^^(1) X (7(23)) X (^x g'^lx^ 



^6 = \/^(''?(23) (7(1) -p- ^(7(23) •a(i))|x 



X7 = y|((?(23)-g-(7(i)-g--^(7(23)-a(i))|x™) 



X8 



3 1 



2 ('^(23) • ^ '?(!) • P + '?(23) • p(7(i) • ^ 



(30) 



- ^p''-g(?(23)-(7(i)j|x™) . (31) 

One can express the \xi) in terms of the Ojlx"*) ^-nd vice versa. Furthermore, we use the most 
general operator structure of the NN t-matrix ttxT'iP^i^ ^ Eg), which is given by 

6 

ttTT'{P,^,Eg) = ^t^^j,,{p,0 ,Eq) U'j((7(2),(7(3),P,p'). (32) 

i=i 

Here the matrix elements t'f^j,i{p,0 , Eq) are scalar functions and the operators Wj((7(2), (7(3),p,^) 
are given in Eq. 
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Next, the 3NF operator V^^^{p,q,p',<f) also allows for an expansion 



ywriP^ 9, P , «') = X] ^u'Ap^ 9, P, q') {^(1) , ^(2) , 0^(3) ,P,q,P,Q) 



,(0 



(33) 



with a strictly finite number of terms. In the context of chiral effective field theory 3NF's have been 
worked out up to N'^LO [jj, llJl, where the expressions can be found. As an example we provide 

nn 

the operators ^li for the 27r-exchange [151 . Il6l | 



'?(2) • Q <?(3) • Q' 



+ Ca(^^ ■ (Q X Q') f(l) • (f(2) X f(3)) 



A + BQ- Q') f(2) • f(3) 



(34) 



where Q and Q' are momentum transfers of nucleons 2 and 3. 



Using 



1^ 



((ti) r|f(2) • f(3)| T) = -6u' 6(-)* I ^ ^ I 

((ti) r|r(l) • (t(2) X f(3))| (t'i) T) = i36 V(2t + l)(2t' + l)(-)*'+l+^ 



i t' T 
t i 1 



1 1 1 

1 1 



2 2 



I, 2 2 



(35) 



(36) 



one obtains 



,(2) 



(1) t^;?^ 



+ %r(p, g,p , ■ ( (p - p ) X (^- g ) 1 , 



(37) 



with Q,Q' expressed in terms of the Jacobi momenta. This defines a subset of the 0; operators 

with the accompaning scalar functions v^'^jj,. 

Now we insert Eqs. (I29|),(l32]) and ([33]) into Eq. ([28]) and obtain 



XI ^tT (P^ 9) Oi (<?(i) , <T(2) , CT(3) , P, I 
1=1 

1 



^ _ £1 _ J_rr2 



+ (-)*4rT'(P' -^(9, 9'), ^g)^«j(o'(2), '?(3)>P, ^))^'23 



2 3 ^ I ^^[HTT'iP^^i^^^)^^Q)'^ji^{-2)^^{3),P,T^{Q,^)) 
m Iml T' j=l 



XI ^tt'T' X ^i'T'i'^'i?^ 9')A^2-P230fc(o'(i), (T(2), <T(3), 7r'(g; q'),q)\x'' 

t' k=l 
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E - ^ - -E-a^ J 



4m'? t' I 

8 

X] "^^t'T (P' Ofc (<?(!) , <?(2) , '?(3) , ^ , 9^ ) I X"") 
fc=l 

^ ~ m ~ T' j=l 

8 

^ </>i'T'(^'' 9")C'fe(a(i), (T(2), <T(3),:p", g'Oix"') 
fc=l 

+ J" / d^qd\"y 

F _ Pi _ _3_„2 / ^ ^ 

m 4m'i t' 



I 

+ (-)*' X] ^«'t(^' 9, -7f(^, 9"). '?{2), 0^(3),^, 9, ^'),q')Pk 

I 

8 

^ Ft,t»T 4't(^'(9"'' 9"), 9"")A'2^l30fc(<?(i), '?{2), '?(3), vr'(g-', g"), 



fc=l 



1 /■ ■ 

m 4m5 T' }=1 



{3),P,P) 



1 



73„' j3„// 



^ m 4m ^ 

+ (-)*' X] ^it'T'(P ' 9, 0^(2), '?(3),P', 9*, -7r(^, q"),q')Pi3 

I 



Pt'f'T' E 4'T'(^'(^"'' r) A'2^'l30fc(a^(l), ^(2), ^(3), ^'(9"', r)!^'"). 



t" k=l 

This can be further simphfied by using 



2 

which turns Eq. (j38p into 

8 

X] "^^tT (P' ((T(i) , CT(2) , CT(3) , P, Ix" 

i=l 

6 



7jOfe(o'(2),<?(3),'?(i)>P,€)(l + '?(23) • o'(i))|x™), 



^ / '^^9' X^Z] (*iTT'(^'^(^'^)'^g)^i('?{2)>^(3),P,7f(^^ 

^ ~ m ~ ITTll T' j=l 



+ -^(^' ^) ' ('^(2) 1 '?{3) ; P> -^(?; 9'))^23) 
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1 

X] ^tt'T' X] '^t'T'(^'(9, 9^), «')(- 2 )C'fc('?(2)> "^(S), o^Ci), 7r'(^, g') 
(l' + a(23)-'a(i))|x'") 



t' I 



K _ El _ ^(,2 

^ m 4m y 
8 

XI ^i'-riP^ ^)Ok{^{i) , CT(2) , (7(3) , :p', g') Ix"" 
fc=i 



1 /" 1 

4-1 J I 



t' " I 

8 

XI '^i'T' (^'' (<7(i) , (T(2) , 0^(3) ,p\q')\x 
k=l 

+ ^ / d^q'd''q"y 



X ^«'t(^' 9' ^(^' ^')^« ('7(1), CT(2), 0^(3), p, g, 7f(g', 9") 

+(-)*' X ^tt'T(^' 9, -T^{q, q"),q")^ii^{i),^{2),^{3),p, q, -T^iq, q"),q")P23 
I 

X Pt'f'T X ^'k^'c?"- 9"), 9"")(-^)Ofc(a^(2), ^(3), ^(1), 7f'(g-', g"), r) 



a + a(23)'- ^(i))lx") 



+ - — ^ — r^X d^p'^^^xAp^'^ ^^q)^0m^^{3),p,p . 

-C/ — — — -e-O^ m/ J 



^ m 4m^ r' " j=l 

d^q'd^q" 



^ m 4m ^ 

X ( X "^"'T' ^(^'j q"),q")^l ('?(!) , '?(2) , 0^(3) , p , q, TT{q, q"),q") 
t' I 

+ (-)*' X'^«'r'(P -^(9''^');9")^«(^(1);'?(2)'?{3)P',^, -7r(g', ^0^23^ 
/ 

8 

X -^i'i'T' X ^t''T' ' 9" ) > ^' ) ( - 2 ) ("^(2) > '^(3) , O'CI) > ' (g' , g ' ) > ^' ) 

(l'+a(23)'^4))lx">- (40) 

Analogous to Eq. dSJ in case of the deuteron we now project from the left with 

(x™'|Ofc(cj(i), C7(2)) '?(3) )P) 9) and sum over m. This leads to the following definitions, which are 
all scalars 

Cki{pq,Pq) = X^^™l^'^('^(i)''^(2)'^(3)'^'^)C'i(o'{i),'?(2),o'(3),p,g)|x"') (41) 

m 

Dkjk'ipq,Pq,f^') = 

11 



m 

Ok' (<?(2) , ^(3) , ^(1) , + ^(23) • a(i) ) Ix™) 

2 ^(^"'lOk (o'(l) , ^(2) , ^(3) , P, Q)Wj (<T(2) , CT(3) , ^, g') 
m 

Ofe'(^(3),'?(2), <?(!), - ^(23) ■ d(i))\xn 

Ekik'{pq,pqp (i,0^) = 

^{x^^Pki^ii), <?(2) , 0^(3) , P, q) 

m 

(^(1) , 5(2) , 5(3) ,P,Q,0,^) Ok' (^(1) , a(2) , a(3) ,0, ^)\x"^) 

Fkjik' im, 0,i> (ii^' ({',1^' 4') = 

XI I ("^(1) ' ^(2) ' ^(3) ' P^ (5(2) ; 0^(3) ,P,f) 

m 

^l{5(l) , 5(2) , 5(3) ,0, Q, 0', ^')Ok' (^(2) , 5{3) , 5(1) , 

(l + a(23)-a(i))|x"^) 

Gkik' {pq, pq0 q,f^') = 
~ 2 ^ ^^'"'^'^ ^"^(1) ' ^(2) ' ^(3) ,p,q)^i (5(1) > 5(2) , ^(3) > p, q, ^) 

m 

Ok'{5(2),5^3),5(i),f,f) (1 + ^(23) • a(i))|x'") 

o'kik' (m, ^V) = 

2^iX^\Ok{5(l), CT(2) , CT(3) , , CT(2) , CT(3) ,p,q,0,^) 

m 

Ok'{S(3).ff(2)5(i)f, q") (1 - a(23) • a(i))|x™) 
^ik,7fc' (m, 0Jqfq",f'q") = 

- 2 ' ^(2) ' ^(3) ' ^^3 ("^(2) > 0^(3) , P, P') 

m 

^Ko'(l) , ^(2) , ^(3) , q, P", ^OOfe' (0^(2) , 0^(3) , , 

(l + a(23).a(i))|x'") 
^fc,7fe' (m, P^.^q^'^'.^"^') = 

2 ^^"^ I (^(1) ' ^(2) ' ^(3) ' ^'^j (^(2) 5 0^(3) , P, p') 
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(l-a(23)-a(i))|x") . (49) 
Using these expressions we end up with the final form of the Faddeev equation 

8 
i=l 

— ^^—^ I d'l' E E 

^ m Am'i T' k'=l 

6 

+ {-)h^r)rj.,{p, -7T{q,q),Eg)D'i,ji,,{pq,p{-7:{q,q)),7f'{q,q)q') 
E Ett'T'(p[/T> (vr'(g, ) 



+ 



E - 




m Am " 




1 


E - 


m 4m " 


E 


1' Sp"d\" 








1 



8 

fc'=i 

m 4my t' 

8 

( E ^«'t(^' 9' ^(^' 9'); E Gkik'iMpqT^iq, q")q', ^^'iq, q")q") 

I k'=l 

8 

+ (-)*' E ^«'t(P' ^, -T^{q, «") E '^'kik'{m,pq{-T^{q , q")W, ^^'i^, q')q" 

I k'=l 
^F,,,T4'^{7f'{q',q"),q") 
t" 

m 4my T' i=i ^ ~ ~ ~ i;;!^ 

8 

E (E'^«'T'(p '^'^(^'^')>^') E ^kjlk'{pq,Pp',pqTriq,q'')q",TT'iq',q")q' 
t' I k'=i 

+ (-)*' E ' 9"-, -v?(9', r), g") 

8 

E ^'kjik'iPQ^PP ^i^Qi-^i^ i^'iQ 

k'=l 
t" 

The first part of Eq. (I50p refers to NN forces only and provides the bulk of the binding energy. 
The remaining parts refer to the action of the 3NF and its interplay with NN forces. 
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It should be stressed, that the momenta only occur in scalar expressions (|4ip -()49p which has to 
be determined only once. For selected examples we refer to the Appendix. 

Due to the expansion of Eq. (j29p the number of coupled equations to solve is strictly finite. 
There at most three 3N isospin states and therefore the number of equations is = 3 x 8 = 24. 
If one neglects charge independence and charge symmetry breaking in the 2N system then N = 8 
(see also p]). Moreover as has been documented in [9, 17 1 only very few out of the 8 components 
are dominant in the expansion of the full '^He state. We expect that a similar reduction takes place 
also for the Faddeev amplitude. 

If one neglects the 3NF only the first part in Eq. (I50p is present. The big question at present, 
however, is to investigate the contribution of the N'^LO 3NF's , which are parameter free and 
therefore the full set has to be treated. 

The final remark refers to symmetry requirements for the unknown scalar amplitudes (p^l^ (p, q) . 
They follow from Eq. (j26p and the expansion given in Eq. (|29|) : 

8 



P23 '^tri-p, Q) = J2 viri-P^ ^ Oi ((?(!), <T(3), (T(2), -p, q)B 



23 IX 



i=l 



= i-YT.^tkp^'1^0,{a^^^,a^2),'^(3),P,q)\xn- (51) 

i=l 

The operators 02,03, O7 and Og are odd under exchange of particles 2 and 3, whereas Oi, O4, O5 
and Oq are even. 

Furthermore, Pislx™') = "Ix™") leading to 



Yl ^tri-P' ((T(i) , a(3) , (7(2) , -p, 



23 IX 



i=l 



E ^^t(p^ 9)0,(^(1), <T(2), a(3),P, q)\xn- (52) 

i=l 

Now we use the symmetry properties of the operators Oj and get 

('ASt^ q)Oi - (k^S i-P^ Q)02 - 49 i-P^ Q)Os + 49 i-p, q)0,+ 

4S{-P, q)05 + 4S i-P^ Q)06 - 4if i-p, q)07 - 4S (-P, ^lOg) Ix'") 

= E0S?(P,91O.|X"^) • (53) 

1=1 
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This is fulfilled if the scalar functions obey 



iT ' 


-p,q) = 


(-)*+v1t(^^'^) 


,rp 1 
Tl 


-P,Q} = 




j.rp 1 
Zl 


-P,Q) = 




tT ^ 


— Tl n) — 
HJ — 


\ ) VtT \P^ 'iJ 


(5), 
tT ^ 


-p,q) = 




(6), 
tT ^ 


-p,q) = 




(7), 
tT ^ 


-p,q) = 


i-U?^{p,q) 


(8), 
tT ^ 


-p,q) = 





(54) 



Our standard manner 



18| to solve the 3N bound state Faddeev equation is by iteration and using 
a Lanczos type algorithm, which requires a choice of initial amplitudes (/'^^^(p, In order to 
guarantee the Fermi character for the 3N system, the symmetry requirements given in Eq. (j54p 
have to be imposed for the initial amplitudes. It is not difficult to explicitely demonstrate that the 
coupled set of equations (jSOp conserves this symmetry requirement during the iteration provided 
the numerics is reliably under control. In Ref. ^ it is documented that the p ■ q dependence of 
the corresponding scalar functions for the full ^He state is rather weak. This suggests to use a 
Legendre expansion for that angular dependence, which would directly guarantee the symmetry 
requirements. 

It might be advantageous to replace the operators Oi by the second form of expansion given in 
Eq. (|29p which is directly related to a partial wave expansion. Ref. [9] documents numerically that 
only very few terms dominate this expansion. It can be easily checked also for the Faddeev ampli- 
tudes which of the two expansions are more efficient, since solutions for realistic forces are available 
in partial wave form and from there the scalar functions (p{p, q) and (j){p, q) can be generated as 
laid out in Ref. 



IV. SUMMARY AND OUTLOOK 

Using an expansion of the Faddeev amplitudes in products of scalar functions cjif^ {p, q) with 
products of scalar expressions in spin and momenta as well as corresponding expansions of the NN 
t-operator and 3N forces, the Faddeev equations for the 3N bound state can be reformulated into 
a strictly finite set of coupled equations for the amplitudes (j)^j^{p,q), which only depend on three 
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variables. The single corresponding equation for three bosons has been solved reliably P, 2]. For 
three nucleons this turns now into a coupled set of equations, which can be solved with modern 
computing resources. 

In case of 3N scattering one standard form of Faddeev equations is [3] 

T|$) = tP\^) +tPGo|$) 

+ (1 + tGo)V^^\l + P)|^>) + (1 + tGo)FW(l + P)GoT|$) . (55) 

We expect that a generalization of the operator expansion for the NN t-operator can be found for 
T|<1>) as well. However, like for the system of three boson, these scalar functions will depend on 5 
variables [sl: \p\, \q\, p-qo, q-Qo, and {% x q) ■ {qo x P), where % is the initial projectile momentum. 
The numerical treatment is thus more demanding with respect to interpolations compared to the 
3N bound state. But again, this has been already mastered P, ^ for 3 boson scattering well into 
the GeV region. 
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APPENDIX A: A SELECTED SCALAR FUNCTIONS 

Below we provide expressions for 2N scalar functions Akk'{p) and Bkjk"{p-,P^)'- 

^ii'(p) = 3 (Al) 
Ai2'{p) = A2v{p) = Q (A2) 
A22'{p) = \p^ (A3) 

Biii{p,p') = 2, (A4) 
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^211 {P, P') = Bi2l {p, p') = 5221 {P, P') = Bu2 {p, P') = Bi22 {p, p') = (A5) 

Bi3iip,p') = {pxp'f (A6) 

Bi4iip,p') = {p + p')' (A7) 

Bi5iip,p') = {p-p'f (A8) 

Biei{p,p') = -2 (pf^-p^) (A9) 

B23iip,p') = -p\py^ - {p- p'f) - \p\p X p'f (AlO) 

B2aMp') = 3(p' +p-p'f - {px p'f - lp\p + p'f (All) 

B2^Mp') = ^{jp" - P- p'f - (P X p'f - \p^ip- p'f (A12) 

^261 = 6# • {p-p')p- {p+p') + 2{pxp'f - ^p\p'^-p'^) (A13) 
Expression for the scalar functions Bij2 are obtained from 52^1 replacing p ^ p' . 

B2i2{p,p') = 3(p • p'f -{pxp'f- ^py (A14) 

B222ip,p') = -8ip- p' {p X p'f (A15) 

B232{p,p') = l{-27{pxp'f+p^p'^{pxp'f + 3p'\px{pxp')f 
y 

+ Sp'^ip' X ip' X p)f + 9ip-p'f{pxp'f) (A16) 

B242(P,P') = ^P'p'\p + P'f-P\p"+P-Pr-P'\p'+P-PT 

+ liP^+p'')iP>^p'f (A17) 

B2,2{p,p') = ^p'p'\p-p'f-p\p"-p-pT-p'\p'-p-pT 

+ lip'^+p'^Wxp'f (A18) 
17 



B2Mp') = 2(p' - p'^kIp'p" - ip- p'f + \{px p'f) (A19) 
The 3N scalar functions Cki(pq,pq) are symmetrical in k and i indices. The diagonal values are: 

Cumm = 2 (A20) 

C22mPQ) = 6 (A21) 

Csampq) = Cumpq) = 2(#x (A22) 

C55{pq,pq) = Cee{pq,pq) = ^^^^'' (A23) 

C77{m,pq) = 2^2^2 (A24) 

CssiplPQ) = 2^'^' (A25) 
The nonvanishing nondiagonal values are: 

C25{pq,Pq) = C2Q{pq,pq) = 2p- q (A26) 

C27{m,pq) = 2^2 (A27) 

C28{pq,pq) = (A28) 

C56 {pq, pq) = C78 {pq, pq) = 2{^-^f (A29) 

C57{pq,pq) = Ce7{pq,pq) = 2p'^p- q (A30) 

C58 {pq, pq) = Ces {pq, pq) = 2q^p -q ( A3 1 ) 

Some other arbitrarily chosen 3N scalar functions Dj^jj^iipq^p'q' ,p"q"), 
D'kjk'{pq,p'q',p"q") and Hkjik'{pq,PP',p'qp"q",P"'q") are: 

D2ii{pq,p'q',p"q") = -3 (A32) 

D222{pq,p'q',p"q") = D243{pq,p'q' ,p"q") = (A33) 

D2u{pq,p'q",p"q") = -3 (A34) 

D'222{pq,p'q',p"q") = (A35) 

H2223{pq,pp',p'qp"q",p"'q") = (A36) 

Hm28{pq,pp',p'qp"q",p"'q") = ip- {p-p') x {p'" - q") 
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- X ((p - p') X p'")) ■ {{p- p') X q") 

- \i0 X {{p-p') X q")) . {{p-p') X p"') (A37) 
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